The necessity of dark matter in MOND within galactic scales 
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To further test Modified Newtonian Dynamics (MOND) on galactic scales - originally proposed 
to explain the rotation curves of disk galaxies without dark matter - we study a sample of six strong 
gravitational lensing early-type galaxies from the CASTLES database. To determine whether dark 
matter is present in these galaxies, we compare the total mass (from lensing) with the stellar mass 
content (from a comparison of photometry and stellar population synthesis) . We find that strong 
gravitational lensing on galactic scales requires a significant amount of dark matter, even within 
MOND. On such scales a 2 eV neutrino cannot explain this excess matter - in contrast with recent 
claims to explain the lensing data of the bullet cluster. The presence of dark matter is detected in 
regions with a higher acceleration than the characteristic MONDian scale of ~ 10~^''m/s^. This is a 
serious challenge to MOND unless the proper treatment of lensing is qualitatively different (possibly 
to be developed within a consistent theory such as TeVeS). 



PACS numbers; 

The standard (ACDM) cosmological paradigm is based 
on Cold Dark Matter (CDM), a cosmological constant 
A, and classical general relativity. Despite its enor- 
mous success and simplicity, competing models have 
been proposed, the main reason being the still unknown 
dark energy component and the undetectability of dark 
matter. To explain the observed fiat rotation curves, 
Milgrom [H proposed MOdified Newtonian Dynamics 
(MOND), based on the relation /(|a|/ao)a = -V<i>N, be- 
tween the acceleration a and the Newtonian gravitational 
field $N- The constant oq ~ 10^^"m/s^ is motivated by 
the acceleration found in the outer regions of the galaxy 
where the rotation curve is flat. When /, assumed to be a 
positive smooth monotonic function, equals unity, usual 
Newtonian dynamics holds, while when it approximately 
equals its argument, the deep MONDian regime applies. 

MOND has been successful in explaining the dynamics 
of disk galaxies; it is less successful for clusters of galaxies. 
It was promoted 01 to a relativistic field theory by intro- 
ducing a TEnsor, a VEctor and a Scalar field (TcVcS). 
TeVeS has been criticised as lacking a fundamental theo- 
retical motivation. Recently, it has been argued Q that 
such a theory can emerge naturally within string models. 

Here we calculate within MOND the defiection angles 
for two generic density profiles and compare them with 
those predicted in standard lensing. We calculate the 
mass of the lenses and estimate the amount of dark mat- 
ter required. We find that despite the alternative gravi- 
tational fall-off, the masses predicted by MOND are very 
similar to those predicted within standard gravitational 
lensing theory. We conclude that MOND within galactic 
scales needs a considerable amount of dark matter. 

We consider a homogeneous and isotropic three-metric 
with the density parameters "tweaked" to the values in a 
MONDian cosmology. The outcome of our lensing anal- 
ysis depends only weakly on the cosmology, for a rea- 
sonable range of cosmological parameters. A different 
background cosmology mainly results in the change of 



the critical surface mass density^]- 

Assuming that the defiection of photons is twice that 
of non-rclativistic particles and that the photon path is 
nearly linear, the deflection angle a as a function of the 
impact parameter b can be written, for a given cumula- 
tive mass profile M{< r), as (see e.g. Q): 



a{h) = - 



ao 



]p2 ^2]3/2 



(1) 



When the function ,f{x) in the integrand is removed, we 
recover the expression of the deflection angle in stan- 
dard lensing. The function f{x) "modulates" this de- 
flection along the path of the particle depending on the 
ratio between the local acceleration, GM{< r)/r'^, and 
the MONDian characteristic acceleration, oq. Wc will 
henceforth use Eq. ^ to calculate the deflection angle. 
In standard lensing f{x) is set to unity, while in MOND 
we first adopt f{x) = x[l + a;^]"^/^. 

Wc compare observations of strong lensing systems 
(which are most often elliptical galaxies) with realistic 
mass profiles. Spherical symmetry is assumed. In addi- 
tion to the "no-dark-matter" interpretation of the rota- 
tion curves in disk galaxies, we assume that in MOND 
the stellar mass content represents the full mass budget; 
the contribution of other baryonic components such as 
gas or dust is minimal in early-type systems. Their char- 
acteristic surface brightness profile can be represented 
by a Hernquist 3-D density profile Q. The cumula- 
tive mass profile is M{< r) = Mr"^ /{r + rh)^, where 
M is the total mass of the galaxy and rh the core length 
scale, related to the projected 2-D half-mass radius by 
Rc = 1.8153 rh. This density model has a logarithmic 
slope (dlogp)/(dlogr) (x — 1 towards the centre, chang- 
ing to —4, as r — > oo. This is our first model. 

The Navarro-Frenk- White (NEW) profile [1] is our sec- 
ond model. The cumulative mass profile diverging loga- 
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FIG. 1: Left: Graphical representation of the lens equation in standard lensing (solid lines) and 
MOND (dashed lines). Each line corresponds to one of the two images of the background source. 
The distant one (number 2) corresponds to the lower set of lines (i.e. a more discrepant result 
between standard and MONDian lensing). The intersection point of the two lines gives the position 
of the source and the total mass (Hernquist profile assumed). Right: NIR HST/NICMOS grey-scale 
image of the lensing system (from the CASTLES database). 



rithmically, we assume a truncation radius rviriai- This 
profile has two free parameters, the core length scale rg, 
and the virial radius. Their ratio is the concentration 
C. Cosmological simulations suggest concentrations 
on galaxy scales to be C ^ 10. Denoting by x the ratio 
X = f/rviriab the cumulative mass function of the NFW 
profile reads 

ln(l+C2;)--^ 



M{< r) = M- 



l+Cx 



ln(l + C) - ^ 



(2) 



i+c 



The lens equation (3 — 6 — a{9)Di^s/ Ds relates the ac- 
tual position of the background source /?, with the posi- 
tion 6 of the images. For a given cosmological model, the 
angular diameter distances from the lens to the source, 
and from the observer to the source, Dls and Dg respec- 
tively, are obtained from the observed redshifts. The 
deflection angle a depends on the mass profile of the sys- 
tem and the impact parameter. A characteristic aspect 
of strong gravitational lensing is that one image appears 
inside the Einstein radius te and the other one outside. 
The difference between MONDian and standard lensing 
lies mostly in the position of the image outside te. 

Figure [T] illustrates our methodology in estimating the 
masses of galaxies from lensing data. HE1104-1805 is 
extracted from the CfA- Arizona Space Telescope Survey 
(CASTLES sample. It consists of a galaxy at red- 
shift zl = 0.73 with a background QSO at zg = 2.32. 
A grey-scale map of the HST/NICMOS F16GW image is 
shown on the right panel, retrieved from the CASTLES 
web-page ^. This is a double system with the image po- 



^ |http://cfa- www.harvard.edu/castles/ 1 



sitions located at 2.09 and 1.10 arcsec on either side of 
the lensing galaxy. The left panel of Fig. [T] shows the 
correlation between the actual position (3 of the QSO, 
and the total mass of the lensing galaxy, assuming a 
Hernquist profile with the projected 2-D half-mass ra- 
dius being equal to the observed half-light radius of the 
lensing galaxy. Each set of lines - dashed (MOND) or 
solid (standard lensing theory) - are the results for each 
image. The compatible solution corresponds to the cross- 
ing of the lines, shown in the figure with a star symbol. 
This gives the true position of the source and the mass 
of the galaxy. For comparison, the values from Refs. [lH 
(for conventional lensing theory) and [l3| (for MOND) 
are given as a shaded region and an arrow, respectively. 

Table I compares our mass estimates with the MON- 
Dian analysis of Ref. [1^ and with the standard non- 
parametric approach of Ref. [Hi (where spherical sym- 
metry is not assumed). The masses are quoted in units 
of IO^^Mq for a ACDM cosmology and, in brackets, for 
the open cosmological model of Ref. [l^l ■ A Chabrier [l3| 
initial mass function is considered for the stellar masses 
quoted from Ref. [lH. The resulting synthetic popula- 
tion, constrained by the photometry of the lensing galaxy 
in the optical (F814W) and NIR (F160W) passbands, is 
used to determine the stellar mass content. The sample 
studied here comprises only double systems to be suitable 
for a 1-D approximation of the lens and serves to show 
the differences between MOND and standard lensing. 

Table I shows a small difference in the mass estimates 
between the two different cosmologies considered here, 
despite their density parameters being quite different. 
This is because the angular distance is mostly unaffected 
by the change in the parameters. There are also some 
noticeable differences between a Hernquist and a NFW 
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FIG. 2; Top: Difference between conventional and MOND masses for a NFW model with C = 10 
(filled dots) and a Hernquist profile (hollow dots). The ratio AM = Mstd — Mmond is shown as 
a function of total (standard) mass (left panel) and the ratio between the average lens separation 
over which lensing masses can be reliably measured, and the observed half-light radius (right panel). 
Bottom: Contribution of dark matter to the total mass budget from a comparison between MON- 
Dian lensing and stellar mass, for a NFW model with C — 10 (filled dots). We also include more 
detailed non-parametric conventional mass estimates of strong lenses from Refs. [iJl and [l7| . 



(C = 10) model for the distribution of mass in the lens- 
ing galaxy. Nevertheless, the differences found are not 
large enough to affect our conclusions. One could always 
argue for a Hernquist profile as this is the model that 
a baryon-only MONDian cosmology would favor, given 
that the projected mass distribution resembles the typi- 
cal de Vaucouleur profile of early-type galaxies. However, 
recent lensing work on clusters, most noticeably the bul- 
let cluster [l5| has been used to postulate a 2 eV neutrino 
which would be important on scales of galaxy clusters, 
not on galactic scales [l6[. We present the NFW profile, 
to illustrate the robustness of our claims in rejecting the 
hypothesis of a 2 eV neutrino. 

The top panels of Fig. [5] compares the mass estimates 
between standard theory and MOND for both density 
profiles: Hernquist (hollow dots) and NFW with C = 10 
(filled dots). The mass differences are shown as a func- 
tion of conventionally calculated mass (left panel) and 
Rlens/R-g (right panel). The difference between the con- 
ventional theory and MONDian predictions stays mostly 



within 10%. This is especially noteworthy in systems 
with Rlens/Re ^ 2. Notice that the lensed images probe 
accelerations slightly above the MONDian threshold. For 
instance, in lens HEl 104-1805 (figure 1), image 2 (right 
panel) is located on the lens plane at a point with a lo- 
cal acceleration of 4.5 x 10~-^°m/s^ (using the MOND 
mass estimate in table 1 for a Hernquist profile), which 
explains why the difference between the solid (standard 
lensing) and the dashed lines (MOND) in the leftmost 
panel is so small. 

The bottom panels of Fig. [2] puts this result in con- 
text with the need for dark matter. The figure com- 
pares MONDian lensing mass with stellar mass (solid 
dots). Our 1-D estimates are compared with more de- 
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tailed non-parametric models from Refs. 
A typical error bar from these estimates is also shown. 
Even though some of the systems can be compatible with 
no dark matter, the MONDian analysis presented here 
finds in most cases the need for dark matter at a level 
around Mdm/-Wstar ^ 0.5-2. Given that the dust and 
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gas content in early-type galaxies corresponds to a frac- 
tion of the stellar mass, we infer the need for dark mat- 
ter even within MOND. Our analysis shows that dark 
matter in early-type systems appears in regions with dif- 
ferent absolute accelerations compared to disk galaxies. 
Hence, a theory with a fixed absolute acceleration (such 
as MOND) cannot explain both early- and late-type sys- 
tems. 

The form of the function f{x), which varies smoothly 
from the deep MONDian to the standard regime is an 
extra source of uncertainty in the MONDian mass esti- 
mates. If f{x) varies too slowly, lingering close to the 
conventional regime for too long, MONDian mass pre- 
dictions are too high, while if f{x) falls quicker to the 
MONDian limit, the need for dark matter would di- 
minish. There is no precise way to determine the ex- 
act form of this function. From galactic rotation curves 
some restrictions can be placed on its form, but there 
still exists a degree of freedom. Varying the form of 
f{x), it was found [l2[ that the predicted masses are not 
affected considerably and that many of the lenses still 
give a high dark matter content. Here, we considered 
two alternatives for the acceleration function, namely 
fix) = x/(l + x) and /(x) = 1 - e^^. The MOND 
mass estimates are lowered by less than 10%. Note that 
one could manufacture a function /(x) such that MOND 
can be successful without dark matter, however such ar- 
tificially made functions would disregard the data from 
rotation curves. 

Another possible source of uncertainty lies in the abso- 
lute value of the acceleration scale ag. One can increase 
ao by a factor 2 and still be compatible with the rota- 
tion curve data [isj . In our case, the mass estimates are 
lowered by about 10%. A combination of a higher qq 
and a shallower function f{x) can result in mass esti- 
mates lower than our fiducial MOND estimates by about 
25% which would still not be large enough to make dark 
matter unnecessary. 

In this paper we have compared mass estimates for 
a set of early-type lensing galaxies using both standard 
lensing analysis and MOND. We used two density pro- 
files, the NFW profile and the Hernquist profile. We 
used the lensing equations to predict the mass of a sys- 
tem from the image positions for a 1-D model (spherical 
symmetry). Besides the standard paradigm ACDM cos- 
mology, other recent alternatives from the literature were 
considered, including the possible solution presented in 
Ref. [3] where the addition of massive neutrinos allows 
a cosmology of (fim, ^^A, i^k) = (0.22,0.78,0) to give an 
acceptable fit to both the CMB angular power spectrum 
as well as the high-redshift supernova data. For our pur- 
poses, any of the cosmologies discussed give very similar 
mass estimates, a result which should not come as a sur- 
prise since the observational constraints mostly impose 
limits on the luminosity and angular diameter scales. 



We tested MOND by looking at a set of strong grav- 
itational lensing early-type galaxies from the CASTLES 
survey. The masses predicted in the framework of con- 
ventional theory are very close to those from MONDian 
lensing, even for galaxies observed out to a few effec- 
tive radii. Comparing the stellar mass content from a 
comparison of the observed optical and NIR photometry 
with stellar population synthesis models we found that 
a very similar amount of dark matter is needed in both 
conventional and MOND analysis. This result is in re- 
markable contrast with the recent attempts to explain 
the lensing data on cluster scales by introducing a 2 eV 
neutrino |16j . This component can cluster on Mpc scales 
but should not cluster on galactic scales to keep the anal- 
ysis of the rotation curves of disk galaxies unchanged. 
However, our lenses, which do require dark matter, are 
studied over length scales comparable to those of the ro- 
tation curve analysis. We therefore conclude that either 
lensing must work in a qualitatively different way within 
MOND (or more correctly the covariant "parent" the- 
ory, such as TeVeS) or dark matter should be considered 
within MOND even on galactic scales. 
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TABLE I: Mass estimates (in W^°Mq units) for ACDM cosmology: (0^, f^A, f7k) = (0.3,0.7,0). The 
masses in brackets correspond to the open cosmology of Ref. (f2m, f^A, fik) — (0.03, 0.36, 0.51). 
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